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Abstract — In this paper, we derive a maximum likelihood 
(ML) decoder of the differential data in a decode-and-forward 
(DF) based cooperative communication system utilizing uncoded 
transmissions. This decoder is applicable to complex-valued 
unitary and non-unitary constellations suitable for differential 
modulation. The ML decoder helps in improving the diversity of 
the DF based differential cooperative system using an erroneous 
relaying node. We also derive a piecewise linear (PL) decoder of 
the differential data transmitted in the DF based cooperative 
system. The proposed PL decoder significantly reduces the 
decoding complexity as compared to the proposed ML decoder 
without any significant degradation in the receiver performance. 
Existing ML and PL decoders of the differentially modulated 
uncoded data in the DF based cooperative communication system 
are only applicable to binary modulated signals like binary phase 
shift keying (BPSK) and binary frequency shift keying (BFSK), 
whereas, the proposed decoders are applicable to complex-valued 
unitary and non-unitary constellations suitable for differential 
modulation under uncoded transmissions. We derive a closed- 
form expression of the uncoded average symbol error rate (SER) 
of the proposed PL decoder with A/-PSK constellation in a 
cooperative communication system with a single relay and one 
source-destination pair. An approximate average SER by ignoring 
higher order noise terms is also derived for this set-up. It 
is analytically shown on the basis of the derived approximate 
SER that the proposed PL decoder provides full diversity of 
second order. In addition, we also derive approximate SER of 
the differential DF system with multiple relays at asymptotically 
high signal-to-noise ratio of the source-relay links. It is shown 
by simulations that the proposed PL decoder in the differential 
DF cooperative system with more than one relay also achieves 
the maximum possible diversity. 



I. Introduction 

Multiple-input multiple-output (MIMO) technology, pro- 
posed approximately a decade ago, revolutionized research in 
the field of wireless communications. By installing multiple 
antennas at the transmitter, benefits like diversity gain and 
spatial multiplexing can be achieved as compared to the single- 
input single-output (SISO) system jTl, Q. It is proposed in 
literature ||3], ||4] that a relaying node can cooperate with a 
source node in order to support the transmissions of the source 
node to the destination. Therefore, the cooperative system 
utilizes distributed antennas to realize a virtual MIMO system 
which can provide benefits of a collocated MIMO system fSl. 

There are two main protocols, namely amplify-and-forward 
(AF) and decode-and-forward (DF), proposed for the coop- 
erative communication systems fS). In the AF protocol, the 
relaying node scales the received data before transmitting it to 
the destination node in order to satisfy a power constraint over 
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the total transmit power The destination requires knowledge 
of the channel coefficients of all links involved in cooperation 
for decoding the data of the source. Whereas, in the DF 
protocol, the relay first decodes the data transmitted by the 
source and then retransmits the decoded data to the destination. 
The destination requires knowledge of the channel gains of 
the source-destination and relay-destination hnks for decoding 
the data of the source. Since the relay cannot decode the data 
perfectly, therefore, erroneous relaying causes significant error 
floor in the performance of the destination receiver This is 
the primary reason why the AF based cooperation has been 
explored in much more detail as compared to the DF based 
system EJ-lIll. 

Differential modulation is useful for the cooperative system 
as it enables the destination to decode the data of the source 
without any channel knowledge of the links involved in the co- 
operation. Therefore, the differential modulation is helpful in 
improving the spectral efficiency of the cooperative communi- 
cation system ifTOl - lfTSl . Differential modulation for a symbol- 
wise DF based cooperative system utilizing uncoded trans- 
missions with one source-destination pair, a single relay, and 
binary phase shift keying (BPSK) constellation is proposed for 
Rayleigh fading channels in ifTOl . lfT2ll and Nakagami-?Ti chan- 
nels in lfT6ll . ifTTl . A maximum-likelihood (ML) decoder for 
differentially modulated binary frequency shift keying (BFSK) 
signal transmitted through multiple orthogonal regenerative 
relays using the symbol-wise DF protocol in the uncoded 
cooperative communication system is found in llT4l . This ML 
decoder considers the possibility of erroneous transmission 
from the relay terminal and maximizes the probability density 
function (p.d.f.) of the received data in the destination terminal. 
In this way, it improves the diversity order of the DF based 
differential cooperative system lfT4l . A low complexity sub- 
optimal piecewise linear (PL) decoder of the binary data is 
also obtained in lfT4l . which performs close to the ML decoder. 
Unfortunately, the ML and PL decoders of the differential DF 
cooperative system obtained in lH4l are not applicable to the 
higher order complex-valued unitary and non-unitary constel- 
lations which are required for increasing the data rate of the 
wireless communication system. In literature ifTSl . differential 
modulation is also used in bidirectional relaying using A/-PSK 
constellation, however, no ML or low complexity PL decoders 
are obtained for this set-up. 

In this paper, we consider a cooperative system with one 
source-destination pair and a single unidirectional relay for 
simplicity. Our main contributions are as follows. 1) We 
derive an ML decoder for the DF based differential coopera- 
tive system utilizing complex-valued unitary and non-unitary 
constellations, and uncoded transmissions. 2) In order to 
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significantly reduce the computational complexity in decoding 
of the differential data, we derive a PL decoder for complex- 
valued unitary and non-unitary constellations. 3) We derive 
a closed-form expression of approximate uncoded symbol 
error rate (SER) of the proposed PL decoder with M-PSK 
constellation. 4) It is analytically proved that the proposed 
PL decoder in the differential DF cooperative system with 
one source-destination pair and a single relay achieves second 
order diversity. 5) An expression of the approximate uncoded 
SER of the differential DF cooperative system with multiple 
relays under asymptotic condition of enw-free source-relay 
links is derived and it is shown by simulation that the proposed 
PL decoder achieves the maximum possible diversity. 

The rest of this paper is organized as follows: In Section II, 
the system model is introduced. Section III derives optimal 
and low-complexity sub-optimal decoders of the differentially 
modulated unitary and non-unitary constellations in the DF 
based uncoded cooperative system with a single relay. In 
Section IV, analytical performance analysis of the DF based 
differential cooperative system with a single relay and M- 
PSK constellation is performed. Differential modulation for 
the DF cooperative system with multiple relays is studied in 
Section [V] The simulation and analytical results are discussed 
in Section VI. Section VII concludes the article. The article 
contains one appendix. 

II. System Model 

Let us consider a cooperative system containing one source, 
one destination, and a single relay as shown in Fig.[T] Each 
node contains one antenna and it can either transmit or receive 
the data at a time. The transmission from the source to the 
destination is performed in two orthogonal phases. In the first 
phase, the source broadcasts uncoded data to the relay and the 
destination. The relay demodulates the data of the source in 
symbol-wise manner and transmits the demodulated symbols 
to the destination in next phase. The source remains silent in 
the second phase in order to maintain orthogonality between 
the transmissions IS). The destination decodes the data of 
the source by utilizing an ML decoder in the second phase. 
We assume a differential cooperative system, where the relay 
and the destination do not require knowledge of the channel 
coefficients of the source-relay, relay-destination, and source- 
destination links for decoding of the data transmitted by the 
source. 

Let in the n-th time interval, the source needs to transmit 
a symbol x[n\ G A, where is a unit-norm Af-PSK constel- 
lation containing the following points: {xi, X2, xs, xm}- 
Before transmission of x[n], the source performs the following 
differential encoding: 

= — l]a;[n], n=l,2,3,..., (1) 
where i;[0] — 1 is an initialization symbol. As |a:[n]|^ = 1, 
therefore, |v[n]|^ — 1. The data received in the destination 
during the first phase in the n-th time interval will be 

ys,d[n] = hs,dv[n] + e^^dW, (2) 
where hs_d represents the circular complex Gaussian channel 
gain of the source-destination link with zero mean and ^ 
variance, and es,d[n] is the signal-independent complex-valued 
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Fig. 1. Cooperative system with a single relay. 

additive white Gaussian noise (AWGN) with zero mean and 
Ns.d variance. During the first phase, the data received in the 
relay in the 72-th time interval can be written as 

VsA^A = hs,rv[n] + es.r[n], (3) 
where hs^r denotes the circular complex Gaussian channel gain 
of the source-relay link with zero mean and a1 variance, and 
es,r[?^] is the AWGN noise with zero mean and N^^r variance. 
It is assumed that the channels hg.d and /ig ,. remain constant 
over at least two consecutive time-intervals n — 1 and n. 
From ([U and (O, we have 

VsA^] = VsA^ - ^x[n] + e'^A^], (4) 

where e', ^[n] = es^,.[n] — esA^ ~ l]a^N ^ CJV {0,2Ns,r) 
is the AWGN noise and CAf (/-t, jy) denotes the complex 
Normal distribution with /i mean and ij variance. It can be 
seen from Q that for given UsA''^ ~ 1] ^^'^ x[n], ys,, M ~ 
CA/" (j/s.r[« — 2iVs r). By maximizing the conditional 

p.d.f. of ?/s,rN given that ys.r[n — 1] and x[n] are known in 
the relay, the ML decoder in the relay can be obtained as 

a;r[n] = arg max Re {y* ^[n]ys An ~ l]x} , (5) 

where Xr[n] G A and = 1. In the second phase, the 

relay differentially encodes Xr [n] into Vr [n] by using ([T]l and 
transmits it to the destination. 

III. Optimal and Suboptimal Decoders in the 
Destination 

In this section, we will derive the ML and PL decoders of 
unitary and non-unitary constellations in the differential DF 
cooperative system. 

A. ML Decoder of M-PSK Data in the Destination 

The data received at the destination from the relay in the 
n-th time interval will be 

yr,dM hr.dVr[n] + CrA^A, (6) 

where hr.d ^ CN {q, ""r d) channel gain of the relay- 

destination Hnk and CrAn] ^ CN{Q,Nr,d) is the AWGN 
noise. It is assumed that hr.d remains constant over at least 
two consecutive time intervals n — 1 and n. From ([T]i and (|6]l, 
we have 

yrAA = yrAn - ll^^rN + e^<j[n], (7) 
where ej,^[n] = er,rf[n] — er.din ~ l]xr[n] is the zero mean 
AWGN noise with 2Nr^d variance. Depending upon the erro- 
neous demodulation of the data of the source in the relay, the 
conditional p.d.f. of yr.d[n] can be written as 



Py.An] (2/|yr,d[n- l],a;N) 

= (1 - e)Py^An] (2/|yr,d[n- l],a;rM = x[n]) 
+ '^Pvr-.dHiylVrA'^ - ^]'^r[n] ^ x[n]), (8) 
where e is the average probability of error 
of the source-relay link. It can be noticed 
from (Ell that py^ ^[n] iy\yr.d[n - '^-], Xr[n] = x[n]) 
CAf (yr,(i[n - l]x[n], 2Nr.d)- From dTji and fill Section III], it 
can be seen that py^ ^[„] (y|?/r.d["- — 1], a^rM 7^ x[n]) denotes 
the p.d.f. of a Gaussian mixture random variable. With this 
observation we have 
PvrAn] {y\yr.dW- l],Xr[n] ^ x[n]) 



27rNr,d{M - 1) 



(9) 



where it is assumed that x[n] ~ Xp, p e {1, 2, Af}, 
Xp G A, is the symbol transmitted by the source. By noticing 
that ^ CAf {ys,d[n - l]x[n],2Ns^d) and owing to 

independence of ys.d[n] and yr.din], an ML decoder can be 
obtained by maximizing the joint conditional p.d.f. of j/s.dM 
and yr,d[n'] as follows: 

x\n] = arg max | ^^Re {yldMysAn - ^]x} 
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+ln ^(1 - g)giv^M<d["l!''-.<'["-i]^} 

N^R<i{y'.dHvr.dln-M^i} 
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It can be seen from ( fTOl i that the proposed decoder requires 
the destination to possess knowledge of the average probability 
of error of the source-relay link e which is a function of the 
average signal-to-noise ratio (SNR) of the source-relay link. 
The relay can estimate the average SNR of the source-relay 
link and feed forward it to the destination. The destination 
can calculate the value of e by using the average SNR of the 
source-relay link and use it for the ML decoding. Since statis- 
tics of a channel vary slowly as compared to the instantaneous 
channel values, the destination requires less frequent updating 
of the average SNR value of the source-relay link. Further, (fTOl i 
is applicable to the A/-PSK constellation, whereas, the existing 
decoders of the differential cooperative system llT2l . llT4l are 
applicable to binary signaling only. For M ~ 2, ( flOl i reduces 
to the existing ML decoder of the binary data lfT2l . lfT4l . 
However, the proposed decoder is computationally complex. 
Therefore, it is desired to have a low complexity decoder of the 
DF based differential cooperative system, which can provide 
maximum possible diversity. 

B. PL Decoder of M-PSK Data in the Destination 

Let ys.d[n] — 0, Vn, i.e., direct link between the source and 
the destination is absent, then a log-likelihood ratio (LLR) 
based symbol-wise decoder can be obtained from (|8) and (|9} 
as follows: 



where p, q 
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+ 92, 
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TABLE I 

Values of T for different values of e used in Fig.O 



e 


10-^ 


10"^ 


10"'' 


10"** 


10""* 


10"" 


T 


±4.9053 


±7.3032 


±9.6148 


±11.9183 


±14.2210 


±16.5236 




Fig. 2. Plots of fit) versus t for different values of e and 16-PSK 
constellation. 
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and 92 = {I- e)e^'^^{''*.-["l?..,d[«-ilx,} ^ 
X gW^Re{,;,M,.,,[„-i].4^ decoder O is 

used as follows to decide about Xp or Xg: ^ ^ 0. Since 
a,ffij52 > 0, it can be shown after some simple algebra that 
^-pi 5: is the same as In ^ ^ 0. 
Hence, we have the following LLR decoder from (firb : 

A^^^ = fit) =ln" 



M-l 



ivb^^ 1 yr.dMyr.d[n - 1] {xp ■ 



(12) 



where t = j^ReAy;j^[n]yr,d[n-l]ixp~Xg)y It 
can be seen from ( 112 b that f{t) clipps to T = 
±ln ((M — 1)(1 — e)/e) for very large and very small values 
of t. Further, it is shown in Table |T] and Fig. |2] for 16-PSK 
constellation that we can approximate f{t) by a piecewise 

linear function as follows: 

( -T, if t < -T, 

fit) « /pL(t) = <^ t, if -T < i < T, (13) 
[ T, if t>T. 

From ( [Tol l and ( fT3T l. we get the following suboptimal PL 
decoder in the destination if the direct link between the source 
and destination is present: 



A, 



to + f?dt). 



(14) 



A, 



where % = j^^^ \y*s,dVAys,dW - 1] {xp ~ .Tg)|. The pro- 
posed PL decoder is applied in a pair-wise manner to the con- 
stellation points for taking a decision of the symbol transmitted 
by the source. For example, for QPSK constellation containing 
the following four signal points: xi,X2,X3,X4, the proposed 
PL decoder decides that Xk,k ~ 1,2,3,4, is the transmitted 
symbol if k^d > Oyi e {l,2,3,4},fc ^ I, where Afe_; is 
calculated from (fl4] i. 

The total number of the real additions and real multipli- 
cations required by the proposed ML decoder dTHt and the 
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Fig. 3. Plots of the total number of real additions and real multiplications 
required by the proposed ML and PL decoders for different A/-PSK constel- 
lations. 

proposed PL decoder ([14) for decoding a symbol belonging to 
the Af-PSK, AI > 2, constellation can be obtained after some 
manipulations as 15M^ + 20M and 33(A/ — 1), respectively. 
We have plotted the decoding complexity (the total number of 
the real additions and real multiplications) of the proposed ML 
and PL decoders for different 7\/-PSK constellations in Fig.O 
It can be seen from Fig. [3] that the decoding complexity of the 
proposed PL decoder increases linearly, whereas, the decoding 
complexity of the proposed ML decoder increases exponen- 
tially with the size of the M-PSK constellation. Moreover, the 
proposed PL decoder significantly reduces the computational 
effort as compared to the proposed ML decoder for M > 4. 

C. ML and PL Decoders for Non-Unitary Constellations 

The non-unitary constellations like Af-QAM are useful for 
increasing the data rate and coding gain as compared to the 
A/-PSK constellation in communication systems. Therefore, 
treatment of the differential modulation based cooperative 
communication system with non-unitary constellation is im- 
perative. For non-unitary constellation like M-QAM, |a;[n]p is 
not necessarily equal to unity, therefore, we need to modify (HJ 
to satisfy the average power constraint as 



v\n — 1] 



x[n\ 



(15) 



where |-| denotes the absolute value. From Q and ( fTSl l. we 
have 



a; 71 



where e' 



x[n~l]\ 
esJn] - esJn - 1] 



(16) 
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|a:[n-l] 

x[n] 



I,(l + - 



a:[n] 
c[n-l] 



As 



CAfiv \n 11 n I 

UJ\ \ys,r[n |a;[„-l]| ' U+ |3.[„_i] 

the relay demodulates x[n] as 



Xr [n] = arg min < In | 1 



therefore. 




where x is an A/-QAM constellation. It can be seen from ([TtI i 
that the differential decoder assumes that x[n — 1] is perfectly 



known to the receiver Since it is not possible for the receiver 
to have perfect knowledge of x[n— 1], the differential decoder 
of ( [TtI i can utilize an estimate of x[n— 1]. 
From (|2]i, (|6]l, and ( flSl l, we can write 



yr,d[n- - 



X n\ 



x[n~ 1]| 
Xr \n\ 



\Xrin - 



where e'^ j^[n] = esM[n] - es,d[n - 1] 
CM (0,{l + \x[n]\y\x[n-l]f)N,., 



11 3:r[Ti] 
' \x,.[n-l]\ 



(18) 



and 



CAf{0,{l + \xr[n]\^/\xr[n~lY)Nr,d)- The joint 
conditional p.d.f. of ys.din] and yr,d[n] can be maximized to 
obtain the following ML decoder of the differential A/-QAM 
data: 



x\n\ = arg max 
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(1 



\x\'- 



l=:r[n-l)l^ 
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E 



|2;,[n-l]r 



1 



M 



i«^.^H-«,,d["-ii| ,^f„'_n ir \ ' 
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(19) 



Since the destination does not have perfect knowledge of x[n— 
1] and Xr[n — 1], it can utilize the estimated values of x[n~ 1] 
and Xr[n — 1] in (fT9] l for decoding the currently transmitted 
data of the source. The estimate of Xr [n — 1] in the destination 
can be obtained by using the following differential decoder: 

2 

yr,d[n - 1] - yrA"^ - 2] I 



Xr[n — 1] = arg min^^ 

xex 



\xr[n-2]\ 



-In 1 




r.d 



(20) 



It will be shown in Subsection IVI-AI through Fig. |6] that use 
of the estimated value of x[n — 1] in the relays and estimated 
values of x[n — 1] and Xr[n — 1] in the destination does not 
lead to error propagation in the performance of the destination 
receiver for different M-QAM constellations. 

After some algebra, the PL decoder of the differential M- 
QAM data can be obtained as follows: 



A, 



uo + fphiur), 



(21) 
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where 



Wo 



\x[n-l]\^ + \x,Y 



\x\n - 



J/s,dN - VsA^ - 1] I 



(l+|.T,|Vk[n-l]|')A^s,d 



= ln 



(l + |xp|Vl^[^-l]|')iV,,d ' 
/ \xr[n - ' ' 



1 + |xqp/|a;r[n - : 

?/r-,dH - 2/r,d[n - 1] 



N.. 



r.d 



\xr[n-l]\ 



(22) 

'l+\Xp\y\Xr[n-l]\'jNr.d 

and fphiur) is given in ( fT3] l. The proposed PL decoder of 
the M-QAM constellation given in d^TT i provides decoding 
complexity which is linear in M. Therefore, the proposed PL 
decoder of the M -QAM constellation also reduces the decod- 
ing complexity as compared to the proposed ML decoder ([19}. 

Remark: We utilize one-way relaying in this paper for sim- 
plicity. The ML and PL decoders obtained in this paper can be 
extended to the two-way relaying based differential DF system. 
However, the extension can be more involved. 

IV. Performance Analysis of the Differential 
Cooperative System with A/-PSK Constellation 

In this section, we will analyze the uncoded symbol error 
rate of the proposed PL decoder of the DF based differential 
cooperative communication system with M -PSK constellation. 

A. Average SER of the Differential DF System with M-PSK 
Constellation 

Let x[n\ ~ Xp, Xp £ A, be the il/-PSK symbol transmitted 
by the source and the destination wrongly decides that x[n] — 
Xq, Xq £ A, is the transmitted M-PSK symbol. The probability 
of error of the destination receiver can be expressed in terms 
of three mutually exclusive events. The conditional uncoded 
pairwise error probability (PEP), given that the channel gains 
of the source-destination and relay-destination links are known 
in the destination, is 

= Pr{to-T < 0\t < -T,x[n]^Xp}Pr{t < ~T\x[n]^Xp} 
+ Pr{to + T < 0\t > T,x[n] = a:p}Pr{t > T\x[n] = Xp} 

+ Pr {to + t < 0, -T < i < T\x[n] = Xp} , (23) 
where Pr{-} represents the probability of an event. Let 7s, d = 
\hs,d\'^ /Ns^d and 7^,^ = \hr,d\'^ /Nr^d be the instantaneous 
SNR of the source-destination and relay-destination link, re- 
spectively. It is shown in Appendix |A] that the conditional 
uncoded PEP of the PL decoder in a differential cooperative 




Fig. 4. Decision regions for A/-PSK constellation. 

system utilizing the DF protocol when hs,d and h^^d are known 
in the destination, will be 

Pe'"'''{hsM,K,d) = Pe^{hs,d) {Pe^iKd) + PeAKd)) 
+ PeAhs,d) {PeAK,d) + Pee{hr,d)) + Peri^sA, hr,d) 
+ Pes{hs,d,hr,d) + Pe,j{hsM,hr,d) + Peia{hsM, K.d) , (24) 

where the terms Pe-{-),i ~ 1,2, ..,6 and P^ {■,■), j = 
7, .., 10 are given in Appendix lAl 

If /is,d and hr.d are the complex-valued Gaussian random 
variables, then js.d and jr,d will be Xi-square distributed with 
the following p.d.f.s 1201 : 

1 L 

P7=,d(7) = — e 

Js,d 

P7.d(7) = ^e"^, (25) 

where ^s,d = <^l.dl^r,d and = crl^dl^rA are the average 
SNRs of the source-destination and relay-destination links, 
respectively. The average uncoded PEP of decoding Xq in place 
of Xp can be obtained as 



E{p-'"-^{hs,d,hrM)} 

Pe l2)P'r.All)P'y,,Al2)d-fl ^72, (26) 



where E {■} denotes the expectation. It is shown in Ap- 
pendix |A] that the closed-form average PEP of decoding Xq 
in place of Xp in the DF based differential cooperative system 
using the M-PSK constellation will be 



Pee) 



P.. 



Pr. 



P.. 



P. 



where the terms Pei , I 



(27) 

1, 2, .., 10 are given in Appendix lAl 
The constellation diagram of an arbitrary A/-PSK constella- 
tion is shown in Fig.lD The decision boundaries corresponding 
to the symbol xi are also shown. From equiprobability of the 
constellation points it can be deduced that 

1 ^ 

Pr [error] = — Pr [errorjxi] = Pr [errorjii] . (28) 
Let Xp — xi — 1 IS Eba transmitted symbol, then the SER of 
the A/-PSK constellation can be approximated by using the 
nearest neighbors approach iBTI Subsection 5.1.5] as 

P^^ P^^'^ ->r P^'-' \ (29) 
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where P^'^'^ and P^"'^ are the average probabilities of error 
of decoding X2 and xm, respectively, in place of the originally 
transmitted symbol xi = 1. We can calculate P^^'^ and P^m^i 
by using ( |27] |. 

B. Diversity Analysis of the Differential DF System with A 
Single Relay and M-PSK Constellation 

Expression of the SER of the PL decoder given in (l27t 
and ( |29] l is very complicated and it is difficult to draw a 
conclusion from it. Therefore, we will derive an approximate 
SER of the proposed PL decoder with i\/-PSK constellation 
by ignoring the higher order noise in this subsection. For 
simplicity, we assume that Xq is the transmitted symbol and 
the destination wrongly decides that Xp was transmitted. It 
can be noticed that because of symmetry of the A/-PSK 
constellation Pe"' " = p^f^-^p. By ignoring the higher or- 
der noise terms in to and t, it can be shown after some 
manipulations that to ^ Af (rc {js^d^qx} j^smI^I"^^ and 

t ^ N ^Re {7r,da;*x} ,7^,^ where x = Xp 

is the symbol transmitted by the relay. After many algebraic 
manipulations, the average pairwise error probability of the 
proposed PL decoder can be obtained by ignoring the higher 
order noise as 

p;--' «/i(7,T) + /2(7,r) + /3(7,r) + /4(7,r), (30) 

where 00 / 

/i(7,r) = (l-e)^ rQ(T—I0] e-^d7,,, 
7s,d Jo V fI 



Xq and Xr 




X 



ls.dlr,d Jj^,j=0 Jy,_a=0 ^/2jV^d\x 

T 1 (™-»gT-r,d)" 

Q 



l«l^7r-,d 



e 

T 

1r,d 

X e "'^■d dwd^r,dd^s,d 



-W - Zq"fs,d 



Is.d 



h{i,T) 



ls.dlr,d J7,,rf=0 Jyr.d=0 VSt^TM I^^I (^'^ ~ 1) 



Q 



-W - Zq^s,d 



M 



1s,d 

e "i-'-^ 



xe dwd-frMdjsM, 
where Q {■) is the Q-function 



(31) 



Eq. (2.3.10)], 



Re{a;<ja;*}, and Zi = Re{xiX*}. It is difficult to solve the 
integration in (|30] | analytically, therefore, the approximate 
average SER of the PL decoder can be obtained numerically 
from ( |29] | and (l30l l. We can use (l30l l for analytically finding 
the diversity of the PL decoder 

Let us assume that the average SNR of all links involved in 
the cooperation is equal and approaching to infinity with the 
same rate, i.e., ^s,d = lr,d = 7s. r = 7 — >■ 00, where 7s^r is the 
average SNR of the source-relay link. From ll23] Section 5.1], 
it can be shown that for M -PSK constellation and large values 
of 7, e oc 4. Therefore, for large values of 7, T oc ±ln7. 

Lemma 1: The following three inequalities are valid for 
integrals of the Q-function: 




d7, (32) 



where Z,ai,bi,ri > 0, and a > 1. 

Lemma [U can be proved with help of the Chernov bound 
Q{x) < e-^'/2,a; > lEH Section 4.L1]. 

Lemma 2: If Z = In 7, 02 7^ 0, 03 is an arbitrary real- 
valued constant, and &2 > 0, then 



and 



1 Ci 

lim — / e e Ta7 « — , 

^■^°°7Jo 7 

1 1 _(z^32)i , C2 

iim — / e ''2T e ^07=—, 

^^■^1 Jo ^/7 7 

e p „z C3 

lim — / e ""STe -1 d'y ~ 

7^007 _/g 7 



(33) 
(34) 
(35) 



if 1/2(2;) = ^jTT/{2x)t 

r(^) 



where ci,C2,C3 are positive constants. 

Lemma |2] can be proved by applying ll24l Eq. (3.471.9)] 
in the left hand side of ( |33]|, (l3j| l. and OS]), then using 
the relations K^{x) « y^Tr/{2x)e~^ ,x >> in ( l33T l, 
)e^-^ in (O, and e « ^ and K^{x) w 

in ( [35] l, where K^{-) is the modified 
Bessel function 1251 Section 9.6.1] and F (•) is the Gamma 
function 1251 Eq. (6.1.1)]. We have also used the fact that ln7 
varies extremely slowly as compared to 7 for 7 — > 00. There- 
fore, for diversity related calculations In 7 can be assumed 
constant relative to 7. 

Let us consider the QPSK constellation for proving the 
diversity of the proposed PL decoder for simplicity. For QPS K 
constellation, Zq = Re {xqX*} is negative, and Zi ~ Re {xiX*} 
can be negative, zero, or positive. Therefore, by using Lem- 
mas [T] and |2] in ( l3Tl i it can be shown that 



7 



lim /i(7,r) (X 



7— >-oo 



C4 



C5 



■ /l (7) , 



lim/2(7,r)cx4 + /2(7), 

7—^-00 'Y'^ 



(36) 



where C4, C5 are constants and /i (7) , /2 (7) are functions 
containing summation terms each decaying at rate higher than 
7~^ for 7 — s- 00. 

By observing the fact that Q{x) is a decaying function of 
X, we have 

,2T(l-e) 



^2tt \x\ Y Jo 

r'oo -j^ 







V7 



e ^7. 



By using Lemmas [T] |2] in ( l37b . it can be shown that 



hm 73(7, r) (X 



£6 

7^ 



/3 (7) , 



(37) 



(38) 



where cg is a positive constant and (7) is a function 
containing summation terms each decaying at rate higher than 
7^^ for 7 — > 00. Similarly, it can be shown that 



lim/4(7,T)cx/4(7) 

7— ^CXD 



(39) 



where /4 (7) is a function containing summation terms each 
decaying at rate higher than 7^^ for 7 — > 00. 

It can be seen from dlB, (ESll, dM), and that 
the average probability of decoding by a PL decoder 

decays as 7"^ at 7 ^ 00 and, therefore, the proposed PL 
decoder achieves a second order diversity. 

V. Differential DF System with Multiple Relays 

Let us consider a general case of N, N > 1, relays 
cooperating along with the direct transmission. It is assumed 
that the relays and the source use + 1 time intervals, in 
order to transmit the data orthogonally to the destination. 

A. ML and PL Decoders for M-QAM Constellation 

It can be shown after some algebra that for the differential 
M-QAM data and N relays the ML decoder will be 



x[n\ = arg max 
xex 



N 



M 



M 



|a;„i [n — 1] I" 



-In 1 



-111' 




(40) 



where is the demodulated symbol in the m-th, m = 

1,...,N, relay in the ri-th time interval, e,„ is the average 
probability of error of the link between the source and the m- 
th relay, is the signal received from the 7Ti-th relay 



by the destination in the ?7-th time interval, and N^^^d is the 
variance of the AWGN noise of the link between the 7Ti-th 
relay and the destination. Since the destination does not have 
perfect knowledge of x[n — 1] and a;,„[n— 1], it can utilize the 
estimated values of x[n— 1] and Xm [n— 1] in ( l40b for decoding 
the currently transmitted data of the source. The estimate of 
Xm [n — 1] in the destination can be obtained by using the 
following differential decoder: 



Xm 77- 



l]=arg mm< 



■2]| 



\x[- 



N, 



-In 1-t-- 



V |x„[7i-2]|7j ■ ^^^^ 

Moreover, the PL decoder for the differential A/-QAM data 
and N relays can be obtained after some algebra as 



N 



A. 



where 



r ni2 



(42) 



Uo — In 
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x[n~l]\ 




(1 + \xg\^ / \x[n - 1 


1') N.,d 



ys.dM - UsA^ ' 



:ln 



(l + \xp\^ /\x[n-lf^N,,d ' 

/ \Xrn[n - + \Xq\^ \ 



(43) 



(44) 



[l + \xp\y\xAn~l]\')Nr^^d 

Tra-i if ^ 

and Tm = ±ln ((A/ - 1)(1 - e„)/e„). 



B. ML and PL Decoders for AI-PSK Constellation 

By substituting \x\ = \x[n]\ = |a:;m[r7]| = 1 in ( |40b and 
after some manipulations, the ML decoder for the M-PSK 
constellation with multiple relays can be obtained as 

x[n] = ai-g max | -^^Re {y* rfM?/s,d[n - l]x} 
xeA I l\s^d 

^ / 1 r 1 

-t- ^ In ( (1 - e„)e^Mj/;™..N!7.„,d["-i]4 



ni—l 



M 



M 



1 ^/ 

i—l.Xi^x 



(45) 
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Similarly, the PL decoder for the M-PSK constellation can be 
obtained from ( |42] | as 



N 



K,q ~t0+Yl /PL(^m): (46) 
m— 1 

Re I y * ^ ,rf[»-]yr„,d[» - 1] (a^p - and /pL(im) is 



where 



given in (144b 

C. Asymptotic SER Analysis of the PL Decoder with Multiple 
Relays and M-PSK Constellation 

Let us assume asymptotically that 7s,r,„ co and ^s,d = 
lr,^.d = 7, where 7s, r„ and 7r„,rf are the average SNRs of the 
links between the source and the m-th relay and between the 
TO-th relay and the destination, respectively. Since 7s, r„, oo, 
and T,n — > oo, i.e., the m-th relay becomes error-free 
under this asymptotic condition. Hence, the asymptotic SER of 
the proposed PL decoder with multiple relays can be obtained 
by assuming that all the relays are error-free. The conditional 
uncoded PEP of the proposed PL decoder ( |46] | with A/-PSK 
constellation and N error-free relays will be 



N 



Pt<^ tp+y^tjn < 0, — oo < < oo|a:[ 



(47) 

By using the results of quadratic forms in complex Gaussian 
variables ll26l in (|47] |. we have 

fe=Ori=0 



X L„ I 7t 



(48) 



where (•) is the generalized Laguerre polynomial ||25] 
pg. 775], /3 = 2Re{x;x}, and 7* = 7s,d + E™=i7r„.d, 
where 7r„,d is the instantaneous SNR of the link between the 
m-th relay and the destination. The distribution of 74 can be 
obtained from (|25]l and lEOl Eq. (2.1.110)] as 



Pit (7) 



-e ^ . 



(49) 



r(iv + 1)7^+1 

By using the series expansion of the generalized Laguerre 
polynomial i24l Eq. (8.970.1)], we can simplify (|48] | as 



a _3\ji\^\ 00 k+N n 



(2|xp-/?)'-(2|xp + /3)\,^ 



A:=0 n=0 i=0 

k+i 



a/V+ji+fc+2i+l ■ ^-^^^ 

The average uncoded PEP of the proposeafL decoder with 
7\/-PSK constellation and N error-free relays can be ob- 
tained by averaging ( fSOl l over 7^. From (|49l l, ( fSOl l, and 1241 
Eq. (3.381.4)], the average PEP will be 

_j 00 k-{-N n 

-EEE"""^"-^ 



r(iv + 1)7^+1 



k=0 n=0 i=a 

(2|x|2 - (2|5|2 +(3yTiN + k + i + 1) 



2Af+n+fc+2i+l)^|j| 



(51) 



^ The average asymptotic approximate 



where c = ^ 
SER of the proposed PL decoder with multiple relays can be 
obtained from ( [29] l and (ISTT i. It can be seen from (ISTt that 
1 1 



Hence, as 7 00, 



Af+1 ■ 



7 



(52) 



(53) 



Therefore, the proposed PL decoder achieves diversity of 
+ 1 in a DF based differential cooperative system with N 
error-free relays. Let iV = 1, then it can be seen from (l53l l 
that the PL decoder with a single error-free relay achieves 
the second order diversity. We have analytically proved in 
Subsection IIV-BI that the PL decoder with a single erroneous 
relay also achieves the second order diversity. Therefore, the 
proposed PL decoder avoids loss in the diversity because of 
an erroneous relaying node. It is shown in Subsection IVI-BI 
and Fig. |9] by simulations that the PL decoder with > 1 
erroneous relays also achieves diversity of + 1. 

VI. Analytical and Simulation Results 

Simulations are performed with Af-PSK and M-QAM 
constellations. The channels of all links are assumed to be 
Rayleigh fading and constant over multiple consecutive time 
intervals. 

A. Performance of the Proposed Decoders for Unitary and 
Non-Unitary Constellations 

It is assumed that % = % ^ = Ir.d = 7. i-c, all links 
involved in cooperation have equal average SNR value. We 
have shown 7 on x-axis in Figs. |5]- 17] Moreover, we have 
considered a DF based uncoded cooperative communication 
system with a single relay and one source-destination pair 
for simulation results shown in Figs. |5} H] In Fig. |5] we 
have plotted the performance of the proposed ML (fTOt and 
PL dll) decoders for differential QPSK, 16-PSK, and 32-PSK 
constellations. It can be seen from Fig. |5] that the ML and 
PL decoders work approximately similar for all constellations 
and SNR values considered in the simulations. We have 
plotted the performance of the proposed ML and PL decoders 
by using the estimated values of x[n — 1] and Xr[n — 1] 
in O, ((111), and (EB for differential 8-QAM, 16-QAM, 32- 
QAM, and 64-QAM constellations in Fig. [6l It can be seen 
from Fig. |6] that the proposed ML and PL decoders also 
work approximately similar for all 7\/-QAM constellations 
considered in simulations. The SER versus SNR performance 
of the proposed ML decoder by assuming that x[n — 1] is 
perfectly known in the relay and x[n — 1] and Xr[n — 1] are 
perfectly known in the destination is also plotted in Fig. |6] 
We call this decoder as 'genie added ML decoder'. It can 
be seen from Fig. |6]that there is no error propagation in the 
performance of the proposed ML and PL decoders due to 
utilization of the estimated values of the previously transmitted 
symbols. 

A cooperative set-up with a single erroneous relay and one 
source-destination pair is considered in Fig. |7]for simulations. 
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Fig. 5. SER versus SNR performance of the proposed ML O and PL □ 

A^nnA^ra fnr HiffWonf AJ-ViK nonctpU urionc 









64-QAM 


1 6-QAM 


8-QAm\ 




X. 32-QAM 











5 10 15 20 25 30 

SNR [dB] 



Fig. 6. SER versus SNR performance of the proposed ML O and PL □ 
decoders utihzing the estimated values of the previously transmitted symbols 
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Fig. 7. Comparison of differential cooperative systems utilizing DF protocol 
with a sub-optimal decoder and proposed PL decoder, and AF protocol 1131 
with 8-PSK constellation. 

In Fig. I2I we have plotted the SER versus SNR performance 
of the differentially modulated 8-PSK constellation in the 



DF based uncoded cooperative system with the proposed 
PL decoder having knowledge of the instantaneous SNR 
of the source-relay link and AF based uncoded cooperative 
system ITSl . The total average transmit power per time interval 
is kept the same in the DF and AF schemes. It can be seen from 
Fig. I2] that the uncoded DF based differential cooperative sys- 
tem with the proposed PL decoder outperforms the same rate 
AF based differential cooperative system utilizing uncoded 
transmissions. For example, a SNR gain of approximately 
1 dB is obtained by the differential DF system as compared 
to the differential AF system at SER=10~^. Full-diversity 
ML and PL decoders of higher order unitary and non-unitary 
constellations in the differential DF system are not available 
in literature. Therefore, we have plotted the performance of 
a sub-optimal decodeiQ that does not have any information 
of the SNR of the source-relay link in Fig. |7] Since the 
sub-optimal decoder does not have knowledge of the SNR 
of the source-relay link, it wrongly assumes that the relay 
is error- free, whereas, the relay actually performs erroneous 
transmissions. It can be seen from Fig. |7]that the sub-optimal 
decoder performs poorer to the proposed PL decoder at all 
SNRs considered in the figure. Moreover, the sub-optimal 
decoder looses diversity due to the erroneous transmissions 
of the relaying node. 

B. Analytical Performance of the Proposed PL Decoder 

We have plotted the analytical approximate average SER 
versus SNR plots of the QPSK, 16-PSK, and 32-PSK constel- 
lations in the differential cooperative communication system 
having one source-destination pair and a single relay in Fig. [8] 
It is assumed that % = % = ^r,d = 7- The approximate 
values of the average SER are calculated in closed-form by 
using dZTl ) and ( |29] l, and by numerically solving the integrals 
in ( |30] | and then using ( [29l ). It can be seen from Fig. [8]that the 
simulation results follow the analytical results satisfactorily at 
all SNR values. Moreover, there is no significant degradation 
in the analysis by ignoring the higher order noise terms for all 
constellations considered in the figure. We have also plotted 
the simulated and analytical performance of the differential 
cooperative system with one relay under the condition that 
7s, d = 7r,ti = 7, 7s, r oo, i.e., the channel between the 
source and relay is error-free for QPSK, 16-PSK, and 32-PSK 
constellations, in Fig. [S] It can be seen from Fig. [8] that the 
proposed analysis closely justifies the simulated behavior of 
the DF based differential cooperative system with error-free 
relay. 

In Fig.|9] we have plotted the analytical approximate asymp- 
totic SER of the differential cooperative system with a single 
source-destination pair, iV = 2, 3 relays, QPSK constellation, 
7s, d = lr,d = 1, and ^s,r 00. The analytical asymptotic 
SER versus SNR values are calculated from (|29t and (ISTT i. We 
have also plotted the simulated SER versus SNR plots of the 
proposed PL decoder for = 2,3, QPSK constellation, and 
7s, d = 7s, r = 7r,d = 7- It Can be noticed that the analytical 
asymptotic SER obtained from (|29] | and ( BTT i is a lower bound 
of the SER of the PL decoder with erroneous relays. From 

'The sub-optimal decoder can be obtained by putting e = in llOt . 
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Fig. 8. Analytical (A SER obtained in closed-fomi by using fZT\ and (29) 
and -k SER obtained by numerically solving the integrals from | |29I and (30), 

and ignoring the higher order noise) and simulated performance of the 

proposed PL decoder with ys,d = 7r,d = 7s,r = 7 and 73^^ = 7,.^^ = 
7i7s,r — > 00, where dB< 7 < 36 dB. An uncoded cooperative system 
with a single relay is considered in the simulations and analysis. 




Fig. 9. Analytical O performance of the proposed PL decoder with 7^ = 
7r d = 7i 7s, r — > 00 and simulated —O— performance of the proposed PL 
decoder with 7s d = 7r d = 7s, r = 7, where dB< 7 < 21 dB and 
TV = 2, 3. 



Fig. |9] it can be seen that for a fixed number of relays, the 
SER versus SNR plot of the PL decoder with erroneous relays 
decays at the same rate as that of the PL decoder with error- 
free relays at high SNR considered in the figure. The diversity 
is defined as slope of the SER versus SNR plot IZTl . Therefore, 
from Fig.|9] it can be noticed that the proposed PL decoder for 
a given number (TV = 2, 3) of error-free and erroneous relays 
achieves the same diversity. We have analytically proved in 
Subsection IV-CI that the diversity (slope of decay of the SER 
versus SNR plot) of the proposed PL decoder with N error- 
free relays is + 1. Therefore, from Fig. |9]and the discussion 
above, it can be noticed that the proposed PL decoder with 
> 1 erroneous relays also achieves diversity of + 1. 
The X-axis in Fig. |9] depicts the average SNR of the source- 
destination link. . 



VII. Conclusions 

We have derived optimal and low complexity decoders 
for differentially modulated complex-valued constellations in 
an uncoded cooperative communication system utilizing the 
DF protocol and multiple relays. Moreover, we have also 
derived expressions of the approximate average symbol error 
rate of the proposed PL decoder It is proved by analysis 
and simulations that the proposed PL decoder achieves the 
maximum possible diversity. 

Appendix A 

Derivation of Uncoded Average PEP with A Single 

Relay 

By considering the erroneous decoding in the relay, we can 
write 

Pr {to-T < Q\t < -T, x[n] = Xp} Pr {t < -T\x[n] = Xp} 
= Pr {to - T < 0|t < -T, x[n] = Xp} 
X ((1 — e)Pr{< < — T|a;[n] = Xp, Xr[n\ = Xp} 
+ePr{t < —T\x[n] = Xp,Xr[n] ^ Xp}) 

= PeAKd) (PeAhr.d) + Pe,{hr.d)) ■ (54) 

Similarly, we have 

Pr{to +T < Q\t > T,x[ri] = Xp}Pr{t > T\x[n] = Xp} 
^Pr{to + T <0\t>T, x[n] = Xp} 
X ((1 — e)Pr{i > r|a::[n] = Xp, Xr[n\ = Xp} 
+ePr{t > r|a;[n] = Xp,a;,-[n] 7^ Xp}) 

= Fe4 {hs4) {Pe, {hr,d) + Pe, {Km)) ■ (55) 

For given hs,d, hr^, x[n], and Xr[n], to and t follow quadratic 
form of Gaussian variates. Therefore, from ll26] Section IV] 
we have Pe^ih^^d) = 1 - g{cp,bp,^s,d)^ Pe2{hr,d) = (1 - 
e)g{bp,Cp,jr,d), Pei{hs4) = g{bp,Cp,"fs,d), and P^^Qir^) = 
(1 - e).g(cp, hp, ■^r,d), where 

1 

r 



00 k k k 
fc=0 n=0 



T{k 



n ■ 



1,2T) 



2"(A; - n)! 
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), (56) 



ai,aj, and 7^ are variables, r(-, •) denotes the incomplete 



Gamma function 1251 Eq. (6.5.3)], i„(-) is the Laguerre 
polynomial 125] pg. 775], bp = 2(2|a;p + XpX + XpX*), and 



2(2|a;|' 



xZ,x 



XpX j • 



Whereas, Pr{t < — T|a;[?7] = Xp,Xr[n\ 7^ Xp} leads to the 
cumulative distribution function (c.d.f.) of the quadratic Gaus- 
sian mixture random variable which can be obtained by 
marginalizing the c.d.f. of the quadratic Gaussian random 
variable |26 Eq. (27)] over Xi,i ^ p as follows: 
Pr{t < — T|a::[n] = Xp,Xr[n] 7^ Xp} 



M 00 



2{M 



-7.,d(2|x|2 



r{k - n 



k=0 n=0 

1,2T) 



4'=fc! 



2"(fc-?i)! 



btjr,d s 



where bi = 2(2|a;p + x*x + XiX*) and Ci 



(57) 

2(2|i|2 - 
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x*x - XiX*). From (|54]i, (|56j and dSTJi, we have Pe^iK^d) = Pi{to+t < 0, -T < t < r|a;[n] = Xp, ^^[n] ^ Xp} 



M 



M- 



M-1 



j'^g{bi,Ci,-fr,d)- Similarly, we can obtain Pes{hr,d) = 

i = \ 
M 

'^^g{ci,bi,jr,d)- Next, we can write 



i = l 



Pi{ta + t <0,-T <t< T\x[n] = Xp} 

= (1 - e)Pi{ta + t< 0, -T <t< T\x[n]=Xp,Xr[n]^Xp} 

+ ePr{to + i < 0,-T < t < T\x[n] = Xp,Xr[n] ^ Xp} . 

(58) 

It can be shown after some algebra that 

Pr{to +i < 0, -T < t < T\x{n\ = a;p,a;r[ri] = Xp} 



(w) / pto{z)dzdw 



Pt\ Xj. [n] — a;p (^) / Pto{z)dzdw 

T 
T 

+ I Pt\xAn]=xp{w) I pto{z)dzdw 



T 
T 



Pt\xr[n]=Xp{w)Fto{-w)dw 



Pt\xr[n]=Xp{w)Fta{-w)dw 

Per {hs,d, hr,d) + —^Peg {hs,d, hr,d), (59) 



1 - e 



1-e 



where Ft„{x) is the c.d.f. of to and Pt\xrln]=xp{x) is the 
p.d.f. of t given that Xr[n] = Xp. By change of variable 
in (|59] l. using 1261 Eqs. (15), (27), (29)], series expansion of 
the incomplete Gamma function 



v—l ^ 

r(«,2/) = (z;-l)!e-^^|j, 

k=0 



(60) 



and then using ll24l Eq. (3.381.1)], we can obtain 

ii .rill J 



and Peg {hs,d-, Ka) = (1 - e) 



CXD k 



A-=0?i=0 



where X>' 



■ii ,mi ,/ 



fc!2''- (fc-n)!4" ' 



mi ,/ 



(mi - 0!7^4"-'=-*i-i7(fc - n + «i + l,4r), 

2 mi!4™i (mi— ;)!2' ' 

oo /c oo mi mi— I 

= E^ E^ . Similarly, we have 

.i,mi, fc=0n=0mi=0 !=0 n=0 



T 
T 



Pt|a;,[n]#2;p (w)^to i-w)dw 



Pt\xr.[n]^Xp {w)Ft„ {-W)dw 



-Peg {hs.d, hr,d) H Peio (^s.rfi ^r^), 

e ■ e 



(61) 



where Pt\xr[n]^xp{x) is the p.d.f. of t given that 2;,.[7i] ^ Xp. 
Since i is the quadratic Gaussian mixture random variable, 
the p.d.f. of t can be obtained by marginalizing Pt\xr[n]=xi{x) 
over Xi,i p as 



a/-i 



-I 2. + ^ + - 



i = i k=On=0 



Pt\xAn]^Xp{v) = < 



■ fe!(fc-n)!2" + '' 
M oo k 



L^{-^), v>0, 



- 2t; + - 



i = i fc=0 n=0 
i^p 



/c!(fc-n)!2"+'= 8 ^ > 

(62) 

From (EB, (|62]i, l26l Eqs. (27),(29)], and El 

Eq. (3.381.1)] we can obtain Peg{hs.dihr.d) = 

M 

tE E ^'Mf'^'m;:!!^*;::^, and (/I,,,, = 



M 



M 



M-1 



EEE^':«'2"-^~'7(-r.+fc+i,2i;, 



i=i fc=On=0 

i^p 



M 



^ ^/ ^ h,n mi,/ zi , 



mi ,^ 



i = l 



, where 7(-,-) is the 



incomplete Gamma function 1251 Eq. (6.5.2)]. 

A close examination of the terms Pe- (•) and Pe (•, •) 
of (l24b reveals that the 2-D integration in (|26t is separable. 
After some algebra and using ll24l Eqs. (3.381.4) 
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